SCALING FEATURES OF SELECTED OBSERVABLES AT RHIC 
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We discuss several observables measured by PHOBOS that show common scaling features in 
Cu+Cu and Au+Au collisions at RHIC energies. In particular, we examine the centrality and 
energy dependence of the charged particle multiplicity, as well as the centrality dependence 
of the elliptic flow at mid-rapidity. The discrepancy between Cu+Cu and Au+Au of the final 
state azimuthal asymmetry (elliptic flow), relative to the initial state geometry of the collision, 
can be resolved by accounting for fluctuations in the description of the initial geometry. 



1 Introduction 



The study of heavy-ion collisions at ultra-relativistic energies allows one to combine experimental 
with theoretical efforts in the understanding of the phase-space properties of strongly-interacting 
matter. At extreme conditions of high temperature and density (compared to normal nuclear 
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Figure 1: Pseudo- rapidity distribution for charged 
hadrons in Cu+Cu (closed symbols) and Au+Au 
collisions (open symbols) at ^/s NN — 200 GeV. 
The Cu+Cu and Au+Au centralities were selected 
to yield similar (Apart) ■ The band indicates the 
systematic uncertainty for Cu+Cu (90% C.L.). 
Errors for Au+Au are not shown. 



TO 



3.0 
2.5 
2.0 
1.5 
1.0 
0.5 



T 



PHOBOS 



t Uf i Li lj * 

* ..•*' Saturation 

200/62.4 

* # * ♦ • ••• 




CuCu AuAu 



preliminary 

PP open symbols closed symbols 



200/J30 



100 200 300 400 

part 



Figure 2: Ratio of mid-rapidity densities as a func- 
tion of -/Vp al -t. Au+Au data are for energies 200/19.6, 
200/62.4, 200/130 (closed symbols) and Cu+Cu data for 
200/62.4 (open symbols). The two solid lines around the 
Cu+Cu ratio indicate additional uncertainty from the es- 
timated trigger efficiency. Ratios for p+p are shown at 
Apart = 2. Model calculations are from Refs. @, @]. 



matter), numerical QCD calculations predict a phase transition to a system dominated by par- 
tonic, rather than hadronic, degrees of freedom. Indeed, one of the important conclusions from 
the discoveries at the Relativistic Heavy Ion Collider (RHIC) is pj that in Au+Au collisions 
an extremely dense, highly interacting system is formed, reaching energy densities much larger 
than ~ 1 GeV/fm 3 , the characteristic scale for the QCD phase transition. 

In this conference proceeding, we will concentrate on a few simple observations extracted 
from data collected by the PHOBOS experiment at RHIC. The observed scaling rules address 
common features in heavy-ion collisions (Cu+Cu, Au+Au) and allow for the comparison with 
simpler systems (d+Au, p+p) in a broad range of collision energies (^/s NN = 19.6 to 200 GeV). 
The presented data are collected by the multiplicity detector covering \rj\ < 5.4 and the near 
mid-rapidity magnetic spectrometer of the PHOBOS experiment (described in detail in Ref. [!]). 



2 Scaling features in heavy-ion data 

We focus on three topics of particle production in Cu+Cu and Au+Au collisions: the scaling 
of the overall charged-particle multiplicity, the factorization of the energy and centrality de- 
pendence of particle production and the connection between the observed final-state azimuthal 
particle distribution and the initial state geometry of the colliding system. 

Before we can address these issues, it is important to quantify how particle production 
depends on the underlying geometry of the colliding species. The spatial overlap of the colliding 
nuclei, determined by the impact parameter, is described as the centrality of the collision. 
Centrality is typically parametrized by the number of participating nucleon pairs (iV part /2), or 
the number of binary nucleon-nucleon collisions (N co u), in the overlap region. Both quantities 
grow with increasing centrality (decreasing impact parameter) of the collision and are calculated 
within a Monte-Carlo simulation of the Glauber model. They are finally related to data via the 
fractional cross-section of the nucleus-nucleus interaction from detailed comparison of measured 
versus simulated charged hadron multiplicity (for details see Ref. 0] and references therein). 



2.1 System-size dependence of particle production 



In fig. [Tj we show the charged hadron dN/dn distributions in Cu+Cu and Au+Au at y/s 



200 GeV. Centrality bins are chosen such that the average number of participants in Cu+Cu is 
approximately equal to that in Au+Au. This comparison, in combination with further studies 
for different centrality selections, leads to a simple scaling rule: If Cu+Cu and Au+Au collisions 
at the same collision energy are selected to have the same (N part ), the resulting charged hadron 
dN/dn distributions are nearly identical, both in the mid-rapidity particle density and the width 
of the distribution. The same is true for dN/dn at 62.4 GeV (not shown, see Ref. 0]). Similar 
findings apply for transverse momentum distributions at both energies [4|. 

2.2 Factorization in energy and centrality 

In studying the centrality and energy dependence we have established a further, more subtle, 
scaling relationship that holds for all of the above mentioned observables. As first described 
in Ref. [dJ for the mid-rapidity density, the increase in particle production per participant with 
increasing N pait is independent of collision energy over the full energy range of RHIC from 19.6 
to 200 GeV. This is illustrated in fig. [2J where we show the ratio of mid-rapidity densities as a 
function of N part relative to the 200 GeV data. All ratios, from 200/130 to 200/19.6, are flat 
within the experimental uncertainty. The factorization is obviously violated in a model, such 
as HIJING pj, that determines overall particle production using a superposition of independent 
"soft" contributions scaled with iV part and "hard" contributions scaled with iV co ii. However, the 
observed factorization is fulfilled in approaches based on the ideas of parton saturation, as e.g. in 
the calculation from Kharzeev et al. for the 200/19.6 ratio (?J. See Ref. 0] for a more complete 
overview on factorization. 

2.3 System-size dependence of elliptic flow 

The elliptic flow at mid-rapidity, V2, as obtained from the angular distribution of particles 
wrt. the reaction plane, is shown in fig. [31(a) as a function of N part for Cu+Cu [3j and Au+Au @] 
collisions at 200 GeV. Its measurement provides important constraints on the hydro dynamical 
evolution of the collision dynamics and gives insight into the connection between initial state 
and final state effects. The initial geometric asymmetry of the collision is typically assumed to 
be given by the eccentricity of the overlap region of the two nuclei. The average eccentricity for 
a certain centrality class is obtained from a Glauber simulation. The eccentricity, £standard; 

of 

the distribution of participating nucleons relative to the reaction plane is commonly calculated 
by fixing the minor axis of the overlap ellipse to be along the impact parameter vector. Since 
the eccentricity for a fixed value of iV part depends on the size of the colliding species, we scale 
out the difference in the initial geometry and compare ^/(^standard) in fig- 0(b) for the two 
systems. The scaled flow in Au+Au is quite flat as a function of centrality showing that the 
flow in Au+Au relatively closely follows the initial state eccentricity. However, comparing the 
scaled Cu+Cu data to the Au+Au data leads to the paradoxical conclusion that for the same 
A^p art , the smaller system (Cu+Cu) is more effective in translating the initial eccentricity into 
the final state anisotropy. A possible resolution to this paradox that produces consistent results 
in Cu+Cu and Au+Au, as shown in fig. [31(c), can be obtained by refining the definition of 
the initial state eccentricity. Our alternative measure of the eccentricity, £ pa rt> is based on the 
observation that ^standard underestimates fluctuations in the actual participants distribution, 
which would generally lead to too small of a mean value. We define £ par t in each centrality bin 
by calculating the eccentricity for each Glauber event relative to the principal axes of the actual 
participant distribution (see Ref. [s|). By construction, e part is positive definite and leads to a 
finite average value even for the most central events. The smaller number of colliding nucleons 
makes the difference between e s tandard an d e P art particularly important for the Cu+Cu system. 
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Figure 3: Unsealed and scaled elliptic flow as a function of TVpart at mid-rapidity (\rj\ < 1) for Au+Au and Cu+Cu 
at ^/s NN = 200 GeV: (a) unsealed t>2, (b) «2/(estandard) 5 ( c ) ^/(epart)- Only statistical errors are shown. 



3 Summary 

We address selected simple scaling observations extracted from Cu+Cu and Au+Au data col- 
lected by the PHOBOS experiment at RHIC. We find that particle production per participant 
nucleon is very similar in Cu+Cu and Au+Au collisions for collision centralities with equivalent 
number of participants. The Au+Au data at RHIC, and to a similar extent also the Cu+Cu 
data, exhibit the factorization of energy and centrality dependence in particle production at 
mid-rapidity. The elliptic flow in Cu+Cu at mid-rapidity is surprisingly large relative to the 
expected average initial state anisotropy, especially when estimated via the "standard eccentric- 
ity". The result can be quantitatively understood, if one accounts for fluctuations in the initial 
collision geometry, which leads to scaling of V2 relative to the "participant eccentricity" in both 
systems. 
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